Dynamic Scaling in One-Dimensional Cluster— Cluster Aggregation 
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We study the dynamic scaling properties of an aggregation model in which particles obey both 
diffusive and driven ballistic dynamics. The diffusion constant and the velocity of a cluster of 
size s follow D(s) ~ s 7 and v(s) ~ s s , respectively. We determine the dynamic exponent and the 
phase diagram for the asymptotic aggregation behavior in one dimension in the presence of mixed 
dynamics. The asymptotic dynamics is dominated by the process that has the largest dynamic 
exponent with a crossover that is located at 8 = 7 — 1. The cluster size distributions scale similarly 
in all cases but the scaling function depends continuously on 7 and 5. For the purely diffusive 
case the scaling function has a transition from exponential to algebraic behavior at small argument 
values as 7 changes sign whereas in the drift dominated case the scaling function decays always 
exponentially. 

PACS numbers: 05.40.-a, 64.60.Cn, 82.20.Mj, 82.70.Dd 
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I. INTRODUCTION 



Both reaction- and diffusion-limited cluster-cluster 
aggregation (DLCA) have been successfully used to un- 
derstand the dynamics of colloidal aggregation |jj . These 
models predict well both the structure of aggregates and 
the growth behavior in dilute particle suspensions as long 
as the dynamics is dominated by Brownian diffusion. As 
the growth of the aggregates proceeds the sedimentation 
of clusters due to gravitation becomes more pronounced 
altering the growth mechanism and cluster structure. 
This was recently observed in experiments [^) . 

The purpose of this paper is to study dynamic scal- 
ing in one-dimensional cluster-cluster aggregation in the 
presence of a competition between diffusion and drift. We 
show that the dynamics at long times is dominated by 
the aggregation process which by itself would lead to the 
fastest growth. The conventional mean-field theory gives 
the correct dynamic exponent for the field dominated 
case but fails when diffusion dominates. The mean- field 
theory also predicts that the scaling function of the clus- 
ter size distribution in the diffusive (driven) case would 
drastically change when 7 (S) changes sign. Such a tran- 
sition is observed for the diffusive case but not for the 
driven one. The dynamic phase diagram shows four dif- 
ferent regions depending on the relative rates of the dif- 
fusion and drift. 

The paper is organized as follows. Section |fj introduces 
the model and describes the algorithm used in simula- 
tions. In Section HI the dynamic scaling is studied us- 
ing mean-field rate equations approach. The mean-field 
results are compared to simulations in section IV. Sec- 
tion |v| concludes the paper with a discussion. 



II. MODEL 



The field-driven cluster-cluster aggregation (FDCA) 
model is defined on a one-dimensional lattice with peri- 
odic boundary conditions, for simplicity. Initially parti- 
cles are distributed randomly on a lattice of L sites up to 
a concentration <p. Sites connected via nearest neighbor 
occupancy are identified as belonging to the same clus- 
ter. The diffusion coefficient of a cluster of size s takes 
the form D(s) = Dis 1 where 7 is the diffusion expo- 
nent and D\ a non-negative constant. The clusters are 
also driven into one direction with a size dependent drift 
velocity v(s) — vis s , which defines the field exponent 6. 

In simulations a cluster is selected randomly and the 
time is incremented by N(t)~ 1 n^ n ] lx , where N(t) is the 
number of clusters at time t and Sl ma x is the maximum 
mobility of any of the clusters in the system at that time. 
The cluster mobility is defined as Q(s) = C v s s + 2Cbs 7 
where C v and Cd are non-negative constants. The choice 
C v = gives normal DLCA. The cluster is moved only 
if x < r2(s)/f2 max , where x is an uniformly distributed 
random number in the interval [0, 1]. The step is taken 
along (against) the field with probability p (g), where 
p = {C v s s + C D si)/n(s) and q = 1 - p. If after the 
move two clusters are in contact, they are irreversibly ag- 
gregated together. Note that time is increased for each 
attempted move. 

Figure [l] shows an example of the dynamics when either 
the diffusion (fig. la) or the drift (fig. 1 b) dominates the 
large time aggregation behavior. The diffusion and field 
exponents are chosen in such a way that at large times 
the largest clusters are the most mobile ones. In the lat- 
ter subfigure, notice the clear breaking of the reflection 
symmetry in the cluster dynamics as the drift begins to 
dominate. Similar behavior is visible in the early-time 
dynamics of the diffusion-dominated case. 
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FIG. 1. An example of the dynamics in FDCA for (f> — 0.1, 
a) 7 = 0.5, S — —1.0 and b) 7 = —1.0, 8 = 0.5. System size 
L — 1000. The time scales are normalized differently. 



III. SCALING ANALYSIS 

Before considering any specific aggregation rules let us 
first represent the well-known mean-field approach. We 
want to compare different dynamical processes in order 
to find the dominating aggregation mechanisms. Denote 
the number of clusters of size s per site at time t by 
n s (t) and the mean cluster size by S(t). The mean-field 
description of irreversible aggregation which neglects spa- 
tial correlations is given by Smoluchowski's equation |^] 



___ = _ y K (i,j) 
at 2 1 ,JJ 
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where the reaction kernel K(i,j) describes the rate at 
which clusters of size i and j aggregate. It is assumed 
to be a homogeneous function K(ai, aj) = a x K(i,j) with 
~ for i <c j. Kernels are classified by /1 [0 : 

/j > (class I), /i = (class II), and /i < (class III). 
Independent of the class the solution scales for mass con- 
serving systems as n s (t) = S(t)~ 2 f(s/ S(t)). In class I 
the aggregation is dominated by the collisions of large 
clusters with large ones whereas the dominant contribu- 
tion in class III comes from the reactions between large 
and small clusters. In class II these two processes are 
equally important. The class III processes can be identi- 
fied from the form of the scaling function since in classes I 



and II fix) ~ x T but in class III f(x) ~ exp(— x as 

Here we concentrate on the scaling function, on the 
polydispersity exponent r, and on the dynamic expo- 
nent z describing the growth of the mean cluster size: 
S(t) ~ t z . The polydispersity exponent in the mean- 
field (mf) is easily found to be r MF = 1 + A in class I. 
Predicting it for class II processes is still a challenge [|. 
However, for all non-gelling systems, i.e. A < 1, the dy- 
namic exponent is related to the homogeneity exponent A 
asz MF = 1/(1- A) §. 

The upper critical dimension, above which the mean- 
field theory is exact, may be calculated once the reaction 
kernel is known j(| . Consider for a moment the aggrega- 
tion of clusters of fractal dimension df in d dimensions. 
For a DLCA kernel Jf D (j,j) ~ (i x / d t -\-j 1 / d f 'j d ~ 2 (P +p) 
(d > 2) the mean-field theory is not exact in any finite 
dimension Jq] but the deviations are negligible already 
in d = 3 pf. In the driven case, if diffusion and ve- 
locity fluctuations are neglected, clusters move ballisti- 
cally. The collision probability of two clusters is pro- 
portional to the product of the mutual cross-section of 
the clusters and the velocity difference between clusters 
K v (i,j) ~ (i 1/d f + j x / d f) d - x \i s _/|. Thus in the mean 
field description the driven system in d dimensions has 
similar scaling properties as the diffusive one in d + 1 
dimensions and therefore the upper critical dimension is 
infinite for both. 

If both diffusion and drift are present the faster dynam- 
ics, as measured by the associated dynamic exponent, 
could be expected to dominate. This is verified by the 
simulation results, discussed in the next section. Thus it 
is adequate to consider the two dynamic processes sepa- 
rately. For example, in one dimension the scaling proper- 
ties of K v necessitate that A = <5 together with [i = 5 for 
S < (class III) and [i = for 6 > (class II). Thus the 
scaling function should drastically change as S changes 
its sign. In one dimension the collision cross-section is 
independent of the clusters' sizes. Thus the above scal- 
ing analysis is directly applicable to the diffusion-limited 
case, too and there should be similar transition between 
the classes III and II at 7 = 0. 

In one dimension the scaling properties of the reaction 
kernels together with z MF = 1/(1 — A) give the mean- 
field dynamic exponent in the diffusive and driven cases 
as z UF = 1/(1 - 7) and z MF = 1/(1 - 5), respectively. 
The strong fluctuations are responsible for the fact that 
the correct exponent is z = 1/(2 — 7) m the diffusive 
case H[§. The dynamic exponent may on the other 
hand be obtained more simply by considering the two 
length scales coming from the two dynamical processes: 
the diffusive length scale Id ~ \FDi and the ballistic 
one l v ~ vt. Naturally the average cluster size is pro- 
portional to the dominant length scale, i.e. S(t) ~ /, 
which together with D(s) ~ s 7 and v(s) ~ s s results 
in z = 1/(2 — 7) an d z — 1/(1 — 8) for the diffusion 
and drift dominated cases, respectively. The simula- 
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tion results presented in section IV confirm these argu- 
ments. Thus the Smoluchowski approach predicts the 
correct dynamic exponent for the driven case even in 
one dimension. If both diffusion and drift are present 
z = max{l/(2 — 7), 1/(1 — 5)} with the crossover at 
5 = 7-1. 

The average cluster size at the crossover can be esti- 
mated by comparing the pairing time (the time required 
for S — > 2S) due to diffusion, t£ , to that due to drift, 
i^ gg . In the diffusive case the pairing time can be ob- 
tained by considering a random walk on coarse-grained 
system with the lattice constant set equal to the aver- 
age cluster radius R |l(]]. In one dimension the cluster 
density on the lattice is p(t) = N(t)/V — (f), where the 
volume V = L/R. A cluster travels a distance of its 
own radius diffusively in time R 2 /D. As it takes on the 
average p~ 2 steps to pair up t^ gg — R 2 /(Dp 2 ). 

For driven clusters the variation in cluster velocities 
is the relevant parameter. Therefore the pairing time 
is of order t* gg = R/(a v p), where <r v = \J (v 2 ) — (v) 2 
is the standard deviation of the cluster velocities. It 
can be calculated from the velocity distribution p(v) — 
sn s \ds(v)/dv\ which gives a v ~ viS s \j ' I2 — I 2 , where 
I a = J dxx Sa+1 f(x) and the approximation comes from 
replacing the sum by an integral. The proportionality 
constant A = \J I2 — I 2 has to be determined numerically 
from simulations since calculating it would require the 
knowledge of the whole scaling function. The crossover 
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FIG. 2. The numerically obtained scaling functions as a 
function of the scaling variable x = s/S(t) for DLCA (7) and 
FDCA (S) at times 10 4 (•••), 10 s (- -), and 9 ■ 10 5 (— ). Sys- 
tem sizes and number of realizations are (5x10, 50), (5x10, 
1000), and (2 x 10 6 , 2000) for 7 = -0.5, 7 = 0.05, and 5 = 0.1, 
respectively. 
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where ro is the elementary particle radius. 

IV. SIMULATIONS 



In simulations the system sizes range from 5 x 10 5 to 
2 x 10 6 , the data are averaged over 50 — 2000 realizations, 
the concentration is usually at <p — 0.1, and random ini- 
tial conditions are used. Neither the initial conditions 
nor the concentration has any effect on the asymptotic 
dynamic scaling properties as was verified by simulations. 
The timescale is fixed by setting Cd = 1 for DLCA and 
C v = 1 for FDCA if not otherwise mentioned. The mean 
cluster size is calculated using both the number (k = 1) 
and weight averages (k = 2) 
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FIG. 3. Dynamic exponent from simulations either as a 
function of the diffusion exponent 7 (*) or the field exponent 
5 (□). The solid line is given by 1/(1 — 5) and the dashed one 
by 1/(2-7). 



s=l 



(3) 



Both averages scale similarly and the number average is 
used in all the figures following. In order to ensure that 
the scaling regime is reached the dynamic exponent is 
calculated using the method of consecutive slopes \ 



We first consider purely diffusive dynamics, i.e., C v = 
0. We obtain an excellent scaling for the cluster size dis- 
tribution using the scaling form n s (t) = S(t)~ 2 f(s/S(t)) 
(Fig. H) an d the known |8||| result for the dynamic ex- 
ponent z = 1/(2 — 7) (Fig. [M. 
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FIG. 4. The scaling functions as a function of the scal- 
ing variable for DLCA for 7 = -0.05, -0.25, -0.50, -0.75 

(from top to bottom) at the times 10 4 (■■ ■), 10° ( ), and 

9 ■ 10° ( — ). System size L — 5 x 10 5 and data are averaged 
over 25 runs except for 7 = —0.05 over 491 runs. Horizontal 
lines show the crossover region exp(l/7) < x < 1 where the 
scaling functions show typical class II behavior. The data for 
various 7- values have been shifted in the vertical direction to 
make the figure clearer. 

The decay of the scaling function near x — depends 
on the sign of 7 and there is a transition from class III 
(7 < 7c) to class II (7 > 7 C ) at 7 C = in accordance 
with the mean-field analysis. However, the transition be- 
tween the algebraic and non-algebraic decay of the scal- 
ing function is plagued by strong crossover effects. This 
is illustrated in figure ^ where the scaling functions are 
presented for several values of the diffusion exponent. 
The crossover behavior is in excellent agreement with 
the mean-field theory according to which the kernels in 
classes I and III show typical class II behavior for inter- 
mediate x values: exp(— l/|/x|) <a;<l||. In our case 
the fi = 7 and the intermediate x region is presented by 
horizontal lines in figure ^. The dynamics for 7 = can 
be solved exactly to establish that the DLCA belongs to 
class II at 7c. The exact result for the cluster size dis- 
tribution is n s (t) — exp(— £)[J s _i(£) — J s +i(£)], where 
T = Wit and I S (T) is the modified Bessel function @. 
This gives f(x) ~ xexp(—Cx 2 ) ~ x (x — > 0), where the 
constant C depends on the average used to calculate the 
mean cluster size. 

As the scaling function decays faster than a power- 
law in class III the polydispersity exponent r is well- 
defined only for 7 > 0. Although the statistics is in- 
sufficient for a direct determination of the relationship 
t(7) the fits to the scaling function show that r increases 
monotonically with increasing 7 so that r = at about 
7 w 0.7. The scaling theory states that for class II 
n s (t) ~ s- T t~ w for 1 < s < S and t — > 00 with the 
scaling relation w = (2 — t)z The exponent w can 
be obtained more accurately from simulations than r. 
A careful analysis of the data shows that w is roughly 



a constant w ~ 1.50 ± 0.05 for 7 £ [0,0.5]. However, 
w cannot be independent of 7 since necessarily w > z 
which diverges when 7 — > 2. By approximating w ~ 1.50 
near 7 = leads to t(^) sw I.5O7 — 1.00, which is zero 
at 70 ~ 0.67 (compare with the actual result above). 
This approximation is consistent with the exact value 
r(0) = -l @. 

Note that the point 70 ~ 0.7 at which the cluster size 
distribution changes from a non-monotonic function to 
a monotonic one is not the same as the transition point 
between the classes 7 C = 0. In the literature it has been 
argued that in two and three dimensions 7 C is negative 
but these arguments rely on the fact the cluster size dis- 
tribution would change to a non-monotonic function at 
the same point p3|. As this is clearly not the case in 
one dimension it is highly probable that j c — in higher 
dimensions, too. 

The corresponding FDCA-simulations are done using 
Cd = 0. Figure H shows also for this case the dynamic ex- 
ponent as a function of the field exponent together with 
the mean-field prediction. The agreement is excellent ex- 
cept for S > 0.3 for which values the asymptotic regime 
has not been reached. 5 = is a special point: all the 
clusters move with the same velocity but the algorithm 
itself causes intrinsic diffusion resulting in the standard 
random walk value z(8 = 0) = 1/2. 

As in the purely diffusive case the cluster size distri- 
bution exhibits scale invariance n s (t) = S(t)~ 2 g(s/S(t)) 
but now with a bell-shaped scaling function g(x) ~ 
exp(-aH"l) as x -> (see figure |). Thus FDCA be- 
longs to class III. No indication of belonging to the class 
II is seen in the range —1.5 < 5 < 0.7 in contradiction 
with the mean-field theory. The absence of the transition 
shows that although the mean-field analysis gives the cor- 
rect dynamic exponent it fails in the case of the scaling 
function. This is not surprising since the spatial fluctu- 
ations expected to be important in low dimensions are 
completely neglected in equation (|l|). Furthermore, for 
5 > the probability for collisions of large clusters with 
large ones is relatively small compared to large-small col- 
lisions since the decisive factor is the velocity difference, 
not the high mobility of large ones. 

The case S = 7 = of FDCA is interestingly enough 
related to a driven diffusive Ising system (DDS). The low 
temperature coarsening in an Ising chain with conserved 
magnetization and subject to a small external force can 
be mapped almost exactly to the diffusion of domains 
with a size-independent diffusion constant flit . The fact 
that the mapping is not quite one-to-one is reflected 
in the behavior of dimers in the DDS. They perform 
long-range hopping which results in another character- 
istic length scale in the problem [jlEfl , As a consequence 
the domain length distribution does not obey the usual 
dynamic scaling for small cluster sizes like it does in the 
FDCA, although the domain size distributions are oth- 
erwise practically the same jl5| . 
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FIG. 5. Average cluster size for various mobilities and con- 
centrations for 7 = —0.5 and 8 — 0.5 in the diffusive Cd = 1, 
Cv = (O), driven C D = 0, C v = 0.05 (□), and driven diffu- 
sive Cd = 1, C v — 0.05 (V) cases. Data are averaged over 50 
runs and system sizes are 10 6 , 5 x 10 5 , and 10 5 for concentra- 
tions 4> = 0.05 (•••), 0.1 (— •), and 0.5 ( — ), respectively. 

Figure ^| shows the crossover from the diffusion domi- 
nated growth to the field dominated one for three differ- 
ent concentrations. Estimating the unknown parameter 
A in equation (§) using the scaling function of the dif- 
fusion limited aggregation for 7 = —0.5 gives A w 0.2. 
Equation (0) gives the crossover sizes 3, 4, and 10 for con- 
centrations <f> — 0.05,0.1, and 0.5, respectively. These 
values agree reasonably well with the simulations as can 
be seen from figure ||. 
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FIG. 6. The phase diagram in one dimension. Roman num- 
bers indicate the class of the aggregation process. Aggrega- 
tion is dominated by the diffusion (light gray), the field (dark 
gray) or by a gelation transition (white). 

V. DISCUSSION 

The results of our study are summarized in figure ^ 
which shows the dynamic phase diagram with four dif- 
ferent regions. The aggregation is dominated by the field 
or the diffusion. At the phase boundary 5 = 7— 1 the two 



processes give the same dynamic exponent. It is unclear 
which one of the aggregation mechanisms determines the 
asymptotic scaling behavior at the boundary. The dif- 
fusive phase is split into two subphases according to the 
dominating aggregation mechanism. The dynamics may 
also be so fast that the systems gels in a finite time. 

Although this paper has considered d = 1 we can also 
discuss the d > 1 case. Here, complications arise because 
the clusters may have a fractal structure. For the field- 
driven aggregation the clusters will in any case become 
anisotropic with a preferred orientation in the field direc- 
tion. We believe both of these complications affect only 
the phase boundaries of the dynamic phase diagram but 
leave its general structure invariant if temporal scaling 
can be assumed. One particular issue is the existence 
of a field-dominated phase with the scaling function be- 
longing to class II. The comparison of the mean-field ap- 
proach and simulations in higher dimensions is left for 
a forthcoming study. The exact location of the phase 
boundaries would be an interesting problem also when it 
comes to applications to experiments. 

In conclusion, we have studied one-dimensional driven 
diffusive cluster-cluster aggregation. We have shown how 
the scaling function depends on the cluster mobilities 
with diffusive or ballistic dynamics, or both. For the field 
dominated case the dynamic exponent can be obtained 
from simple mean-field calculations, which together with 
the simulation results may be used to obtain the phase 
boundaries in the dynamic phase diagram. This shows 
four different phases in the aggregation depending on the 
relative strength of the diffusion and the field. 
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